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THEORETICAL  LOADING  AT  SUTERSONIC  SPEEDS  OF  FLAT  SWEPT-BACK 
WINGS  WITH  INTERACTING  TRAILING  AND  LEADING  EDGES 
By  Doris  Cohen 


SUMMARY 

By  the  method  of  superposition  of  conical  flows,  the  load  distribu¬ 
tion  is  calculated  for  regions  of  a  long,  rectilinear,  swept— back  wing 
behind  the  points  at  which  the  Mach  lines  from  the  trailing— edge  apex 
intersect  the  leading  edge.  It  is  found  that  a  good  approximation  to 
the  load  can  be  obtained  by  the  application  of  a  fairly  simple  correc¬ 
tion  factor  to  the  two-dimensional  subsonic  distribution.  The  similarity 
to  two-dimensional  flow  is  used  to  derive  expressions  for  the  loss  of 
lift  behind  the  Mach  lines  from  the  tips  of  the  wing. 


INTRODUCTION 

In  a  previous  paper  (reference  l),  formulas  were  derived  by  the 
superposition  of  lineari*zed  conical  flows  (reference  2)  for  the  loading 
of  a  thin,  flat,  swept— back  wing  of  rectilinear  plan  form  traveling  at 
supersonic  speeds,  under  the  limitation  that  the  Mach  lines  from  the 
trailing-edge  apex  did  not  cut  the  leading  edge.  Ahead  of  these  Mach 
lines  the  flow  is  known  to  follow  the  conical  pattern  identified  with  the 
triangular  wing  having  the  same  apex  angle.  By  superposing  an  appropriate 
distribution  of  other  conical  flows  behind  the  trailing  edge  to  cancel  the 
lift  in  the  wake,  the  triangular-wing  loading  can  be  corrected  within  the 
region  behind  the  trailing-edge  Mach  lines  to  conform  to  the  Kutta  condi¬ 
tion  for  a  subsonic  trailing  edge.  This  method  is  described  more  fully 
in  reference  1. 

When  the  region  of  influence  of  the  trailing  edge  includes  a  part  of 
the  leading  edge  and  the  stream  ahead  of  it,  the  mutual  interaction  of  the 
two  regions  results  in  a  flow  which  is  susceptible  to  treatment  by  the 
conical-flows  method  only  at  the  expenditure  of  considerable  effort.  The 
flow  in  the  tip  regions  is  particularly  difficult  to  calculate.  A  method 
of  cancellation  of  external  lift  by  supersonic  doublet  distributions,  pro¬ 
posed  in  reference  3,  offers  some  promise  of  providing  a  solution  in  the 
tip  region  if  the  pressure  at  all  points  inboard  of  the  tip  Mach  cone  is  ^ 
known  from  other  sources.  The  method  appears,  however,  to  be  in  error  as 
applied  in  the  neighborhood  of  the  leading  edge  or  of  a  raked— out  tip. 
Formulas  of  considerable  simplicity  have  been  derived  by  Heaslet,  Lomax, 


2 


NACA  ™  1991 


and  Spreiter  (reference  4),  but  only  for  a  particular  family  of  slightly 
tapered  wings  with  a  filleted  trailing  edge  and  plan  forms  very  slender 
relative  to  the  Mach  cone. 

Since  no  satisfactory  method  appeared  to  be  available  for  investi¬ 
gating  the  general  problem  of  the  rectilinear  swept— back  wing  with 
interacting  leading  and  trailing  edges  —  that  is,  wings  of  high  sweep 
or  aspect  ratio,  or  wings  of  moderate  sweep  and  aspect  ratio  at  low 
supersonic  Mach  numbers  —  the  work  of  reference  1  was  extended,  with 
some  simplifying  assumptions  for  the  tip  regions,  to  cover  such  cases. 

The  principal  object  of  the  investigation  was  to  obtain,  if  possible,  an 
indication  of  the  manner  in  which  the  flow  in  the  outboard  regions  of  a 
long  narrow  wing  approaches  the  subsonic  two-dimensional  type  of  flow 
indicated  by  simple  sweep  theory. 

The  calculations  showed  that  the  section  loading  assumes  the  general 
shape  of  the  subsonic  two-dimensional  loading  at  the  very  station  at 
which  the  trailing-edge  Mach  line  intersects  the  leading  edge.  Further 
investigation  showed  that  everywhere  behind  the  Mach  lines  from  the  point 
of  intersection  and  farther  outboard  the  load  could  be  closely  approxi¬ 
mated  by  a  fairly  simple  adaptation  of  the  two-dimensional  formulas. 

The.  following  material  outlines  the  method  of  obtaining  the  load 
distribution  in  the  region  of  leading—  and  trailing— edge  interaction  by 
the  superposition  of  conical  flows  and  presents  the  formulas  for  approxi¬ 
mating  the  results  by  a  single  step  based  on  correction  of  the  two- 
dimensional  formulas.  A  corresponding  approximate  correction  to  the  load 
near  the  tips  is  also  derived. 

CALCULATION  OF  LOADING  BY  SUPERPOSITION  OF  CONICAL  FLOWS 

Procedure 

Figure  1  shows  the  central  portion  of  a  long,  narrow,  swept-back 
wing  divided  into  various  regions  by  the  trailing-edge  Mach  lines,  the 
Mach  line  arising  at  the  intersection  of  the  trailing-edge  Mach  line 
with  the  leading  edge,  and  so  on.  In  the  superposition  process  to  be 
described,  these  Mach  lines  mark  the  regions  of  influence  of  the  flow 
fields  used  in  successive  cancellations  of  lifting  pressure  outside  the 
wing  boundaries. 

In  region  I,  the  streamwise  component  u  of  the  uppe2>-siarface  per¬ 
turbation  velocity  —  to  which,  in  linear  theory,  the  loading  is  propor¬ 
tional  —  is  given,  as  previously  mentioned,  by  the  known  solution  for  the 
triangular  wing  swept  behind  the  Mach  cone 


u^(a) 


^  Up 
</  m^-e^ 


(1) 
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where  uq^  the  constant  velocity  along  the  center  line  a=0,  is 


mVa 

Uo  - - T=“ 

P  E(  Vl-m2) 

(See  Appendix  A  for  explanation  of  the  symbols,) 


(2) 


The  flow  over  the  remainder  of  the  wing,  however,  is  affected  by 
the  reduction  of  the  lift  to  zero  behind  the  trailing  edge .  The  primary 
induction  effects  on  the  wing  are  given  in  reference  1  in  the  form  of 
two  corrections  to  the  basic  perturbation  velocity:  a  symmetrical 
correction^ 


F( 
°  K( 


(3) 


With 


cp  = 


and  the  integral  with  respect  to  a  of 

(Au)  =  -  iia  cos“i  (l-a)(ta-fflt)  -  ( m^.-a ) ( l-t„  ) 
«  (l-mt)(ta-a) 


(4) 


where  u^  =  (du^/da)da. 

These  induced  velocities  appear  on  the  wing  surface  (in  regions  II, 
III,  •  •  . )  as  corrections  to  the  basic  velocity  u^.  Beyond  the  actual 
surface  of  the  wing  they  represent  a  discontinuity  in  pressiire  which 
cannot  be  supported  in  the  free  stream  and  must  therefore  be  canceled  at 
all  points  ahead  of  the  leading  edge.  Cancellation  of  these  pressure 
differences  will,  in  turn,  induce  additional  corrections  to  the  perturba— 
tion  velocities  in  region  III  and  all  portions  of  the  wing  aft  of  that 
region.  Again,  pressure  differences  will  be  introduced  in  the  wake  region 
and  new  corrections  will  have  to  be  made  which  will  affect  the  loading  in 
region  IV  and  all  more  rearward  portions  of  the  wing.  For  the  present  no 
tip  location  will  be  specified,  except  that  it  is,  of  course,  outboard 
of  the  intersection  of  the  first  trailing-edge  Mach  line  with  the  leading 
edge.  The  tip  region  will  be  considered  after  the  result  of  the  trailing^ 
edge  “  leading--edge  interference  effect  is  determined, 

^Equation  (3)  as  originally  prescribed  in  reference  1  (equation  (51))  is 
incorrect.  An  errata  sheet  has  been  issued  to  correct  this  error. 
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The  initial  distribution  of  lift  to  be  canceled  ahead  of  the  wing  (in 
region  II  extended)  is  the  result  of  superposing  an  infinite-number  of 
conical  flow  fields  originating  at  the  various  points  on  the  trailing  edge. 
To  cancel  this  distribution  of  lift  by  the  superposition  of  conical  flows 
requires  that  each  of  the  lift  distributions  corresponding  to  a  single 
value  of  a  in  equation  as  well  as  the  lift  due  to  (Au)q  (equa¬ 

tion  (3))^  t)e  canceled  separately  by  an  infinite  number  of  conical  fields. 

The  single  conical  field  of  (Au)q  will  be  considered  first.  The 
velocity  field  to  be  superimposed  ahead  of  the  leading  edge  to  cancel  the 
velocities  (Au)o  induced  in  the  plane  of  the  wing  by  the  symmetrical 
solution  (equation  (3))  can  be  built  up,  as  shown  in  the  sketch  on  the  left. 


of  overlapping  constant-velocity  sectors  having  one  edge  along  the  leading 
edge  of  the  wing  and  one  along  the  extended  ray  t©  from  the  apex  of  the 
trailing  edge.  The  magnitude  of  the  constant  velocity  on  each  element  is 
[d(Au)Q/dtQ]  dtQ  or,  from  equation  (3)> 


^o  <^^o 

K(  yi-^2 


(5) 


Similarly,  for  each  oblique  wake  or  trailing-edge  element  a,  (see 
right-^and  sketch),  a  canceling  field  can  be  built  up  ahead  of  the  leading 
edge  by  the  superposition  of  sectors  bounded  by  the  leading  edge  and  by 
rays  ta  from  the  apex  Xa,ya  element  a  and  having  a  constant 

velocity  of  magnitude 
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d(Au) 


a 


N  ‘i'ta  - - - 

“  «  ‘‘at. 


COS 


-1 

(l-]24.)(tg^-a) 


dt, 


(from  equation  (4)). 

Elementary  Solution  for  the  Region  Ahead  of  the  Leading  Edge 


In  general,  the  elementary  solution  required  for  the  cancellation 
of  pressure  in  the  plane  of  the  ving  ahead  of  the  leading  edge  is  one 
that: 


1.  Provides  constant  streamwise  velocity  over  an  infinite  sector 
hounded  on  one  side  hy  the  leading  edge  of  the  wing  (extended)  and  on 
the  other  hy  an  arbitrary  ray  extending  outward  into  the  stream  from 
some  point  X]^,y5  on  the  leading  edge.  (See  fig.  2.) 

2.  Induces  no  vertical  velocity,  or  downwash,  on  the  wing. 

3.  Induces  no  lift  except  on  the  wing  and  within  the  sector 
described  in  condition  (l). 

At  first  glance  these  conditions  would  appear  to  he  satisfied  hy  the 
oblique  solutions  used  at  the  trailing  edge  in  the  previous  work,  if 
properly  oriented  with  respect  to  the  wing,  and  one  might  expect  the  same 
form  of  solution  to  apply.  In  reference  1,  however,  it  was  pointed  out 
that  the  downwash  connected  with  this  solution  remained  constant  over  the 
wing,  only  provided  that  the  wing  area  did  not  include  the  line 
y  =  constant  extending  downstream  from  the  apex  of  the  element.  In  the 
case  of  the  leading-^dge  element  this  condition  is  violated  (fig.  2)  and 
an  additional  term  is  needed  to  bring  the  downwash  to  zero  throughout  the 
area  of  the  wing  affected  hy  the  element. 

The  solution  applicable  to  this  case  has  been  given  by  Lagerstrom 
(reference  2).  The  u  component  of  the  velocity  in  the  plane  of  the  wing 
is  as  follows: 


u 


r.p. 


’  cos-1 


2m 

( l+m)ta. 


>J  (tg^-m)(l+t^) 


(6) 


In  equation  (6),  u-^  is  the  constant  streanwise  perturbation 

velocity  over  the  element,  and  t-^  is  |3  times  the  inclination,  with 
respect  to  the  stream,  of  the  ray  from  the  apex  x-j^,yi2  of  the  element 
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through  the  point  x,j  at  which  the  pressvire  is  being  computed.  Thus, 


■tb  = 


3(y-yb) 


When  the  correction  is  being  made  for  the  symmetrical  trailing-edge  ele¬ 
ment  only,  ta  is  replaced  by  tg  and  u-t,  =  [d(Au)o/d.to]  dtg* 

For  brevity,  the  two  parts  of  the  correction  function  will  be 
referred  to  as 


C(t„)  =  r.p.  cos 


(l+m)(t-b-ta) 


E(t„)  =  r.p. 


( l+m)t, 


(t»-ffl)(l+t„) 


These  functions  and  the  induced  velocity  (equation  (6))  are  plotted 
against  t-jj  in  figure  2  for  typical  values  of  the  other  parameters, 


Leading-Edge  Corrections  to  the  Local  Lift 

Applying  equation  (6)  to  the  cancellation  of  the  symmetrical-wake— 
correction  velocities  (Au)o  ahead  of  the  wing  resiilts  in  the  following 
induced— velocity  increment  at  any  point  x,y  on  the  wing: 


(Au^) 


^  d^  +E(to)]  dto 


vrhere  Tq  is  that  value  of  tQ  for  which  t-j^  =  ”1^  designates  the 

most  rearward  leading-^dge  element  containing  the  point  x^y  within  its 
Mach  cone.  In  terms  of  x  and  y. 


-  =  m(x+Py) 

°  (x+(3y)-(l+m)cQ 


d(Au)o 


Integration  of  -tt — -  C(tQ)dtQ  is  not  possible  by  elementary  means, 
dto 

For  values  of  m  equal  to  or  greater  than  0.4,  a  series  expansion  can  be 
used  for  (Au)q  (Appendix  B)  and  a  satisfactory  expression  found  for  the 
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first  term  of  the  integral  (equation  (lO)).  Calculations  for  m  =  0.2 
indicate  that  integration  by  series  may  not  be  siifficiently  accurate 
when  m  is  small.  In  such  a  case,  the  first  part  of  the  integration 
may  be  performed  graphically. 


The  second  part  of  equation  (10)  can  be  integrated  in  closed  form  as 
follows : 


”v7  dto  °  °  jtmt(l+m)K’(mt)  v  mx-Py 

■^o 

where  K  and  E  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind,  respectively,  of  the  modulus 


KE(k,\|f)-^rp’(k,t) 


(12) 


2mt(l-Tn) 

(l-mt)(TQ+mt) 


(13) 


and  F(k,\|r)  and  E(k,>|f)  are  the  corresponding  incomplete  integrals,  with 
the  argument 


^  =  sin  ^ 


V  2to 


(1^) 


The  corresponding  procedure  for  each  oblique  trailing-edge  element 
a  results  in  an  expression  for  [d(Au)'b/da]  da  similar  to  equation  (lO). 
Again,  the  first  term  cannot  be  integrated  in  closed  formj  the  second  term 
becomes 


/^S(Au)a 

i  -^a 


K(ta)  dta 


“4m  Ug^V  1-a 

Jt2(i+in)a 


/■ 


( mj.-a )  ( x+Py )-( l+m)m(.c, 
mx-py 


X 


(l+a)(T„-a) 
( l-a)(mj--a) 


KE(ka,>lra)-«(ka,^lra) 


where 


”4; 


KE(k^,ilrQ)^(k^,^r^) 


(15) 


m(  mj.-a)  ( x+Py  )-e,m|.CQ(  1+m) 
(  mj.-a )  ( x+Py  )-m^CQ(  1+m) 


(16) 
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In  region  III  (fig.  l),  the  total  leading-edge  correction  to  the  per— 
turbation  velocity  at  any  point  x,y  is 

(Au)b  .  (AUb)„  .  r*"  s*  (20) 

"Where  aQ*  corresponds  to  the  most  rear'ward  trailing—edge  element  the 
velocity  field  of  which^  reflected  in  the  leading  edge,  can  affect  the 
point  x,y. 

In  region  IV,  the  cancellation  of  the  velocity  increment  (Au)^  at 
the  trailing  edge  must  be  taken  into  accoimt.  To  do  this  in  any  rigorous 
way  is  obviously  impractical.  Ho"wever,  the  velocities  to  be  canceled  are 
small,  it  is  possible  to  estimate  the  effect  of  their  cancellation  in 
many  cases.  Beyond  region  IV*  the  velocity  distribution  cannot  be  deter¬ 
mined  in  practice,  except  very  approximately,  by  the  conical— flows  method. 

Calculated  Results 

Load  distributions  have  been  calculated  by  the  foregoing  method  for 
three  combinations  of  taper,  sweep,  and  Mach  number  as  follows: 

Untapered  Tapered 

m  =  0.2  0.4  0.4 

114.  =  0.2  0.4  0.6 

These  values  of  m  and  represent,  by  virtue  of  the  Prandtl-Glauert 

transformation,  a  variety  of  practical  sweep  angles  at  Mach  numbers 
between  1  and  2;  as  for  example,  0.2  would  be  the  value  of  m  for  a  wing 
with  63°  sweep  of  the  leading  edge  at  a  Mach  number  of  1.075^  or  75^  sweep 
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at  a  Mach  numher  of  I.25.  Similarly,  m  =  0.4  would  correspond  to  45^ 
of  sweep  at  M  =  I.08,  60^  at  M  =  1.22,  or  75°  at  M  =  I.80.  The 
trailing-edge  sweep  angles  at  these  latter  Mach  numbers,  if  mt  =  0.6, 
are  34°13S  49°>  and  68^,  respectively. 

In  figure  3,  the  m  =  0.4  wings  are  shown  at  M  =a/2,  so  that 
A  =  cof"!  0.4.  Keeping  the  same  angle  of  sweep,  the  m  =  0.2  case  is 
depicted  (fig.  4)  as  representing  the  same  wing  at  a  lower  Mach  number 
(P  =  1/2).  Figures  3  and  4  show  also  the  spanwise  stations  at  which  the 
load  distributions  were  calculated. 

The  resiilts  of  the  calculations  are  presented  (figs.  5>  6,  and  7)  in 
the  form  of  values  of  P  (Ap/qa)  where  Ap  is  the  local  lift,  q  the 
free— stream  dynamic  pressure,  and  a  the  angle  of  attack.  In  linear 
theory  this  quantity  is  connected  with  the  streamwise  component  of  the 
upper-surface  perturbation  velocity  u  by  the  relation 


^  ^  k  — 

qa  Va 

The  various  components  of  the  lift,  as  calculated  by  the  superposition 
method,  are  shown  individually. 

Figure  5  presents  the  resulting  lift  distributions  at  the  two 
stations  of  the  tapered  wing  shown  in  figure  3«  Section  A— A  contains  the 
intersection  of  the  trailing-edge  Mach  line  with  the  leading  edge,  so  that 
the  leading-edge  correction  at  the  leading  edge  is  zero.  At  points 
farther  back  along  the  leading  edge,  as  at  Py/co  0.8,  the  correction 
is  minus  infinity.  However,  it  is  seen  to  increase  to  a  small  positive 
value  within  a  fraction  of  the  chord  length  at  Py/co  =  0.8. 

At  both  stations  it  is  necessary  to  estimate  the  effect  of  cancella¬ 
tion  of  the  leading-edge  correction  at  the  trailing  edge  to  satisfy  the 
Kutta  condition.  Cancellation  would  be  effected  by  means  of  oblique  ele¬ 
ments  of  the  type  used  previously  (equation  (4))  in  canceling  lift  at  the 
trailing  edge,  so  that  the  induced  lift  on  the  wing  has  the  same  general 
shape  as  the  primary  oblique  trailing-edge  correction.  It  therefore  may 
be  presumed  to  fall  along  a  modified  inverse-cosine  curve  to  zero  at  the 
boundary  of  the  region  affected.  Locating  this  boundary,  the  twice- 
reflected  Mach  line  from  the  trailing-edge  apex,  we  are  able  to  draw  in  a 
satisfactory  estimate  (dotted  cirrve)  of  the  correction  to  bring  the 
pressure  once  more  to  zero  at  the  trailing  edge. 

The  untapered  wing  with  the  same  sweep  (m  =  0.4)  relative  to  the 
Mach  lines  is  shown  in  figure  6  with  the  load  distributions  calculated  at 
the  same  stations. 

• 

Foiir  section  lift  distributions  are  presented  (fig.  7)  for  m  =  0.2. 
At  3y/co  =  0.15  only  the  rear  60  percent  is  influenced  by  the  subsonic 
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trailing  edge.  The  reflection  of  this  influence  in  the  leading  edge 
alters  the  pressure  over  the  rear  40  percent  of  the  section.  At  section 
B-B  the  leading-  and  trailing-edge  interaction  affects  the  entire 
section.  A  further  reflection  of  this  effect  in  the  trailing  edge  must 
he  estimated. 

At  section  C— C  the  reflection  of  the  leading-edge  correction  in  the 
trailing  edge  covers  the  whole  extent  of  the  chord  and  any  estimate  of 
its  magnitude  would  he  necessarily  arbitrary.  Also,  a  second  pair  of 
reflections  must  he  taken  into  account.  The  final  pressure  distribution 
has  therefore  been  drawn  as  a  hand  within  which  the  true  curve  may  he 
shown  to  lie.  Its  height  is  the  error  introduced  at  the  trailing  edge 
by  the  first  leading-edge  correction,  except  very  near  the  leading  edge, 
where  an  infinite  negative  correction  is  known  to  he  introduced  by  the 
second  leading-edge  correction.  The  calculations  were  also  carried  out 
f*or  Py/cQ  =  0.45.  The  margin  of  uncertainty  was  found  not  to  have 
increased  by  any  appreciable  amount.  (See  fig.  7(d.).) 

It  may  be  seen  from  the  figures  that,  except  in  the  vicinity  of  the 
leading  edge,  the  magnitude  of  the  leading-edge  correction  is  small.  The 
part  due  to  reflecting  the  oblique  trailing-edge  correction  is  roughly 
the  same  fraction  of  the  total  leading-edge  correction  as  the  oblique 
trailing— edge  correction  is,  at  the  leading  edge,  of  the  total  trailing- 
edge  correction.  This  fraction  is  less  than  10  percent  in  the  examples 
calculated  and  the  resulting  contribution  to  the  leading-edge  correction 
is  negligible.  The  labor  involved  in  computing  the  loading  can  in  such  cases 
be  considerably  reduced  by  omitting  the  calculation  of  the  effect  of 
reflecting  the  oblique  trailing-edge  correction  in  the  leading  edge. 


APPLICATION  OF  TWO-DIMENSIONAL  FOEMOLAS 
TO  CALCULATION  OF  LOAD  DISTRIBUTION 


Load  on  Midspan  Portion  of  Wing 

It  is  apparent  from  the  calculated  results  that,  whenever  the  plan 
form  and  the  Mach  number  are  such  that  the  trailing-edge  Mach  line  inter¬ 
sects  the  leading  edge,  the  load  distribution  behind  the  Mach  lines  from 
the  point  of  intersection  (regions  III  and  beyond,  fig.  l)  resembles  in 
shape  the  well-known  incompressible  load  distribution  over  an  infinitely 
long  flat  plate.  However,  the  quantitative  agreement  does  not  appear 
good,  particularly  in  the  case  of  the  tapered  wing,  until  the  load  distri¬ 
butions  over  sections  normal  to  the  stream  are  examined,  when  a  simple 
proportionality  of  the  curves  is  observed.  In  order  to  determine  the 
factor  of  proportionality  it  is  noted  that,  if  the  curves  are  similar, 
each  can  be  characterized  by  the  strength  of  the  singularity  at  the 
^leading  edge.  An  explicit  expression  for  this  quantity  can  be  obtained 
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for  regions  III  and  of  the  swept-back  wing^  and  the  ratio  of  this 
quantity  to  the  corresponding  quantity  in  the  two-dimensional  flow  is  the 
required  correction  factor  to  bring  the  two  curves  into  agreement. 

The  leading-edge  singularity  in  the  loading  on  the  swept-back  wing 
is  initially  that  in  the  triangular— wing  loading-  Introduction  of  the 
leading-edge  corrections  to  the  load  in  region  III  reduces  the  strength 
of  the  singularity  there  through  the  terms  R(to)  and  R(ta).  (The  inverse 
cosine  function  is  always  finite.) 

It  is  convenient  to  define  the  strength  of  the  leading-edge  singular¬ 
ity  as  the  coefficient  of  (mx— py)^/^  in  the  velocity  distribution.  In 
the  region  ahead  of  the  trailing-edge  Mach  line,  this  coefficient  is, 
from  equation  (l). 


reducing  to 


at  the  leading  edge. 

From  equations  (12)  and  (15)>  decrements  to  this  coefficient  may  be 
derived  for  the  portion  of  the  leading  edge  behind  the  intersection  with 
the  trailing— edge  Mach  line,  as  follows: 


and^  for  each  value  of  a  from  0  to  that  value  a©  which  makes  Tg^ 
eqiial  to  one. 


dAC 

-  da  = 

da 


-4m 


'(l+a)(Tg^-a) 

(l-a)(mt-a) 


KE(kg,Mro)  -EF(kg,to)  |  (23) 

This  expression  will  be  shown  later  to  hold  approximately  for  some 
distance  outboard  of  its  theoretical  limit  of  applicability,  that  is, 
beyond  region  IV  (fig.  l). 
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For  convenience^  we  may  define  the  nondimensional  coefficient 

In  figure  8,  the  quantity  Pa(x)  is  plotted  against  x/c©  for  the 
wings  previously  discussed.  Additional,  terms  enter  into  the  coefficient 
at  the  values  of  x/cq  indicated  hy  cross-markings  on  the  curves,  where 
successive  reflections  of  the  Mach  lines  off  the  trailing  edge  intersect 
the  leading  edge.  Evaluation  of  these  terms  has  not  been  attempted; 
however,  the  asymptotes  of  the  curves  for  the  untapered  wings  are  known 
from  simple  sweep  theory,  and  it  is  apparent  from  the  figure  that  the 
magnitude  of  the  omitted  decrements,  within  the  practical  range  of  aspect 
ratios,  must  he  small.  The  curve  for  the  tapered  wing,  without  the  addi¬ 
tional  corrections,  goes  to  zero  where  the  wing  tapers  to  a  point,  at 
X  =  3co*  However,  the  complete  curve,  with  the  effect  of  the  successive 
reflections  included,  would  have  a  finite  value  at  that  point.  It  can 
be  shown  that  the  successive  corrections  enter  with  zero  derivatives  of 
successively  higher  order  so  that  the  error  caused  by  neglecting  one  ,^is 
initially  small.  The  curve  in  figure  8  is  believed  valid  to  the  extent 
of  the  solid  line. 
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The  subsonic  perttirbation  velocity  on  an  infinitely  long  flat  plate 
moving  normal  to  its  surface  with  circulation  r  is 


_r_ 

JtC 


(26) 


where  tj  is  the  distance  to  the 
leading  edge,  expressed  as  a 
fraction  of  the  chord  c.  If 
the  section  of  the  plate  coin¬ 
cides  with  a  section  of  the  swept— 
back  wing  taken  perpendicular  to 
the  stream  (x  constant),  as  in  the 
sketch,  then  the  chord  length  is 


I  [mx-^nt  (x-Cq)] 

and 


mx-Py 

mx-m^Cx-CQ) 


(27) 


(28) 


Substituting  for  t]  in  equa¬ 
tion  (26),  we  find  that  the 
coefficient  of  (mx— Py)*"  is 


which  reduces  at  the  leading  edge  to 


(29) 


7TC 


mx-m|.(x— Cq) 


Equating  this  coefficient  to  the  coefficient  for  the  swept— back  wing  gives 


X 

TtC 


_ 1 _ 

V  mx-mt(x-CQ) 


Ca  +  (AC)o 


(30) 


Then  the  two-dimensional  approximation  to  the  load  coefficient  on  the 
outboard  portions  of  a  swept— back  wing  is  given  by 


l4 
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or 


u 

Va 


a(x) 


[py-^Xx-^„)]c„ 

[  mx-^ii,- (  x-Cq  )  ]  ( mx-Py ) 


(31) 


^  =  4a(x)  /  rP7-g4:(x-^o)3Co 

qa  J  [inx-mt(x-Co)](mxH3y) 


(32) 


The  closeness  with  which  the  foregoing  procedure  predicts  the  theo¬ 
retical  loading  over  swept— hack  wings  is  indicated  hy  figures  S,  10,  and 
11,  where  the  previously  calculated  load  distributions  are  con^ared  with 
those  calculated  by  equation  (32).  Even  in  the  case  of  the  highly 
tapered  wing,  the  agreement  is  seen  to  be  good. 

At  the  most  inboard  section  of  the  m  =  0.2  wing  (fig.  7(q*))  there 
is,  of  course,  no  agreement  over  the  forward  portion  where  the  flow  is 
essentially  conical,  but  behind  the  reflection  of  the  trailing-edge  Mach 
line  in  the  leading  edge  (at  60— percent  chord)  the  adjusted  two- 
dimensional  theory  appears  to  hold  fairly  well.  At  the  next  section, 

B“B,  the  agreement  is  very  good.  At  sections  C-C  and  D— D,  the  curve 
calculated  by  the  two-dimensional  theory  lies  within  the  band  prescribed 
for  the  loading  by  the  conical— flow  calculations.  Since  the  discrepancy 
between  the  two-dimensional  loading  and  the  correct  theoretical  distribu¬ 
tion  is  already,  at  section  B^,  less  than  the  width  of  the  bands  in 
figures  7(c)  and  (d),  and  must  diminish  to  zero  at  infinity,  it  may  be 
supposed  that  the  connected  two-dimensional  curve  is  at  least  as  satis¬ 
factory  an  approximation  to  the  correct  curves  at  sections  C— C  and 
D— D  as  at  section  B-B,  and  more  satisfactory  than  can  be  obtained  by 
the  conical-^lows  method. 

The  load  distributions  derived  by  single  sweep  theory  are  included 
(figs.  9(b),  10(b),  and  11(d))  to  show  the  magnitude  of  the  plan— form 
effect  and  also,  in  the  case  of  the  untapered  wings,  the  curves  which  the 
load  distributions  are  approaching  with  increasing  distance  from  the 
plane  of  symmetry. 

In  figures  10(b)  and  11(b),  comparison  is  also  made  with  results 
of  the  slende3>^ing  theory  of  reference  4.  In  the  second  case 
(fig.  11(b)),  the  Mach  line  from  the  trailing  edge  of  the  root  chord,  in 
the  modified  plan  form  of  reference  4,  crosses  the  section  at  12.6 
percent  of  the  chord,  causing  the  cusp  which  appears  in  the  load  distri¬ 
bution  calculated  for  that  wing.  It  is  interesting  to  note  that  the 
loadings  calculated  by  the  two  methods  agree  fairly  well  even  when  the 
wing  violates  the  slender  plan— form  assumption  (m=0.4). 
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Calculation  of  Tip  Effect 

In  the  conical-flows  method,  the  effect  on  the  load  distribution  of 
the  presence  of  the  wing  tip  is  calculated  by  canceling,  outboard  of  the 
tip,  each  of  the  various  components  of  the  lifting  pressure  previously 
calcinated.  To  do  this  exactly  when  the  trailing-edge  Mach  line  and 
leading  edge  intersect  involves  a  prohibitive  amount  of  work.  Approxi¬ 
mate  procedure  s  have  been  developed  (Appendix  C),  but  the  calculations 
remain  tedious,  while  the  accuracy  is  in  some  respects  imsatisfactory. 

A  simpler  approximation,  based  on  the  assumption  of  cylindrical  flow 
beyond  the  tip  location,  was  therefore  attempted.  The  formulas  are 
developed  for  the  tips  alined  with  the  stream,  but  the  fundamental 
procedure  may  be  applied  to  raked  tips  as  well. 

It  is  assumed  that  the  velocity  distribution  to  be  canceled  in  the 
p^n®  of  the  wing  outboard  of  the  tip  is  an  extension  of  the  velocity 
distribution  calculated  for  the  tip  section  along  lines  parallel  to  the 
leading  edge.  For  this  purpose  the  approximate  load  distribution  given 
by  equation  (32)  is  used,  still  further  sinqjlified  by  assuming  a  to 
remain  constant  at  its  value  at  the  leading  edge  of  the  tip  section.  The 
error  involved  in  this  assuii5)tion  may  be  estimated  by  referring  to 
figi^e  8.  The  maximum  range  of  abscissa  to  be  used  is  1.0;  the  change  in 
ordinate  within  that  range  —  where  the  influence  of  the  wake  is  effective  — 
is  less  than  10  percent.  (Where  the  wing  is  tapering  to  a  point  and  o 
is  changing  more  rapidly,  the  tip  region  is  so  small  that  the  entire  cal¬ 
culation  of  tip  effects  would  probably  be  omitted.) 

The  assuiii)tion  of  constant  0  results  in  a  failure  to  cancel  exactly 
the  lift  along  the  tip.  The  assvunption  of  cylindrical  flow,  while  reason¬ 
able  for  the  untapered  wing  (compare  figs.  10(a)  and  (b),  for  example) 
would  appear  to  be  too  drastic  for  the  tapered  wing,  where  neither  the 
chord  nor  the  loading  remains  constant.  However,  it  has  been  shown  in 
reference  1  that  the  major  part  of  the  tip  effect  results  from  the  cancel¬ 
lation  of  the  infinite  pressure  along  the  leading  edge,  and  this  part  will 
be  accurately  calculated.  The  effect  of  the  residual  lift  on  the  rearward 
portion  of  the  tip  section  and  in  the  stream  will  be  shown,  by  numerical 
examples,  to  be  small. 

The  distribution  of  perturbation  velocity  at  the  tip  station  y  =  s 
with  the  simplification  of  constant  0  is,  from  equation  (31)^  approxi¬ 
mately 


u(xc,s) 


=  0, 


[ps-^(xc-Co)]  Cq 

[mxc-mt(xc-xo)  J  (mxc-Ps) 


(33) 
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where  Xc,s  are  the  coordinates 
of  a  point  on  the  tip  and  Oq  is 
the  value  of  a  at  x  =  3s/m. 

This  expression  may  be  more 
conveniently  written  in  terms  of 
the  variable 


which  is  the  ratio  of  the  distance 
of  xc^s  back  of  the  leading  edge 
(see  sketch)  to  the  tip  chord  ct* 
Since 


equation  (33)  ^7  written 


(36) 


If  the  velocity  distribution  jf 

u  is  assumed  to  be  constant  beyond 
y  =  s  along  lines  parallel  to  the 

leading  edge,  it  may  be  canceled  by  the  superposition  of  conical  constant- 
velocity  elements,  having  one  edge  along  the  tip  and  the  other  parallel  to 
the  leading  edge,  with  apexes  displaced  along  the  tip  by  increments  in 
6g.  The  velocity  induced  at  a  point  x,y  by  each  such  element  would  be 
(reference  2) 


jt 


m+tg+2mtg 

tg-m 


(37) 


where 


P(y-s) 

X-Xc 


(38) 


and  Ujj  is  the  velocity  on  each  sector. 

Following  the  procedure  used  in  deriving  equation  (36)  in  reference  1, 
the  corresponding  equation  may  be  written  for  the  velocities  induced  by 
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canceling  the  assumed  cylindrical  loading  : 

Au 
Va 

where  x^^s  is  the  intersection  of  the  Mach  forecone  from  x^j  with  the 
“tip* 

If  the  distances  of  x,j  and  Xq^s  hack  of  the  leading  edge^  meas¬ 
ured  as  fractions  of  the  tip  chord,  are 

and 

("o-t)  (‘‘D 

it  can  he  shown  that 

(l+m)(x-XQ)  =  (42) 

from  which  equation  (39)  can  he  written  (with  the  substitution  for 
u(xc,s)  from  equation  (36)) 


where 


X  =  ^ 


the  taper  ratio 


In  integrating  equation  (43),  three  cases  must  he  distin^ishedi 
(1)  5  <1  (always  true  for  the  imtapered  wing),  (2)  1  <  I  <5^^  (when  the 

point  x,y  lies  behind  the  trailing  edge  of  the  fictitious  untapered  wipg 

through  the  tip  chord),  and  (3)  ^  (o*  possible  condition  for  some 


points  near  the  trailing  edge  of  a  highly  swept  or  tapered  wing). 
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In  the  first  case  (£  <l) 


where  K  and  E  are  the  complete  elliptic  integrals  of  modulus 


and  E(k*,ti)  Q-re  the  incomplete  integrals,  with  the  com¬ 

plementary  modulus 


and  argument 
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In  the  second  case  fl<  I  < 


mt-my 


Au 

;r-(x,y) 

Va 


[Ml(k,\|f2)-EB’(k,t2)] 


(44b) 


^x,y)  =  -  -  [KE(k•,^^3)-{K-E)F(k^^|r3)]  (44c) 

%NmK  /m.-m  \ 

Vi  (^1-0 


where 


Alox}^  the  Mach  line  from  the  leading-odge  tip,  all  three  equations 
reduce  to  the  value 


Au-^  = 

Vo.  ./iSCl" 


(45) 
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By  the  procedure  just  described,  approximate  cancellation  of  all 
pressure  differences  outboard  of  the  tip  has  been  effected,  but  the 
pressures  induced  by  such  cancellation  now  violate  the  condition  of  zero 
lift  in  the  wake.  Cancellation  of  the  induced  pressure  differences  in 
the  wake  can  be  accomplished  in  an  approximate  manner,  as  in  reference  1, 
by  making  use  of  the  known  value  of  the  induced  pressure  at  the  trailing 
edge  of  the  wing,  but  assuming  the  entire  error  to  originate  at  the  tip 
leading  edge . 


On  the  trailing  edge  of  an 
untapered  wing,  |  =  1  and 

"N  ^  fefK  (1^6) 

itv  mX 

There  is  no  corresponding  simplifi¬ 
cation  for  the  tapered  wing. 

The  constant— velocity  elements 
to  be  superposed  in  cancellation  of 
this  tip— induced  velocity  will  be 
defined  by  rays  of  slope  tin/p 
from  the  leading-edge  tip.  (See 
sketch. )  Their  apexes  will  be  at 
the  intersection  xt^yt  rays 

with  the  trailing  edge,  and  they 
will  be  bounded  by  the  rays  tm  on 
one  side  and  the  trailing  edge  on 
the  other.  The  initial  sector, 
bounded  by  the  tip  Mach  line  tm  =  — 1 
will  have  the  constant  velocity  Au* 
given  by  equation  (45).  The  total 
induced  correction  to  the  velocity  at 
any  point  x,y  will  be 


A2u(x,y)  =  - 


Au*  cos' 


(l-iiit)(t*+l) 


( 1-tj^)  ( t-^ )-( mf^m)  ( l-'t) 


(^7) 


where,  if  x*,y*  are  the  coordinates  of  the  intersection  of  the  tip  Mach 
line  with  the  trailing  edge. 
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t*  =  (48) 

and 

t  .£<Z2i1  (1.9) 

X 

Kie  integral  in  equation  (^7)  is  most  conveniently  evaluated  by  plotting 
the  inverse  cosine  against  Au,  obtained  from  equation  (46)  for  the 
untapered  wing,  or  equations  (44b)  or  (44c)  for  the  tapered  wing,  for  the 
points  xt,yt  on  the  trailing  edge.  The  two  functions  are  related 
through  tjn^  which  is  defined  by 


P(yt-s) 


(50) 


It  is  difficult  to  assess  the  accuracy  of  the  formulas  proposed  for 
calculating  the  tip  losses  in  the  loading  since  no  correct  values  are 
available  for  comparison.  In  figure  12  two  streamwise  section  load  dis¬ 
tributions  calculated  by  equations  (kh)  and  (47)  are  compared  with  those 
obtained  by  the  simplified  conical— flows  methods  described  in  Appendix  C. 

The  latter  calculations  are  correct  alon^  the  tip  Mach  lines,  but  are 
known  to  fall  below  the  correct  theoretical  value  at  more  rearward  points. 
The  tip-induced  velocity  to  be  canceled  at  the  trailing  edge  cannot  be 
calculated  with  significant  accuracy  by  the  approximate  method  employed, 
so  that  no  estimate  could  be  made  of  the  load  curve  within  the  region 
affected  by  such  cancellation.  The  cylindrical— flow  approximation  probably 
represents  as  good  an  estimate  as  can  be  obtained  at  present  of  the  theo¬ 
retical  load  distribution  in  that  region. 


CONCnroiNG  EEMARES 

Beyond  the  spanwise  station  at  which  the  Mach  lines  from  the  trailing 
edge  apex  intersect  the  leading  edge  of  a  swept— back  wing,  the  section 
loading  takes  on  the  general  form  of  the  two-dimensional  subsonic  load  dis¬ 
tribution.  A  good  approximation  to  the  loading  in  this  region  can.  be 
obtained  by  applying  an  easily  calculated  correction  factor  to  the  two- 
dimensional  flat— plate  distribution. 

The  two-dimensional  approximation  appears  to  be  in  good  agreement 
with  the  theoretically  more  accurate  loading  calculated  by  the  conical- 
flows  method  for  tapered,  as  well  as  untapered,  wings.  In  the  case  of 
tapered  wings,  the  proposed  method  is  able  to  predict  the  increased 
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loading  in  the  outhoard  regions  (similar  to  the  effect  found  in  subsonic 
flow) . 

The  same  concept  of  utilizing  two-dimensional  loading  can  he  employed 
to  obtain  an  estimate  of  the  loss  of  loading  in  the  zone  of  influence  of 
the  wing  tip^  when  the  tip  lies  wholly  within  the  trailing-edge  Mach  cone. 
This  sinplification  is  especially  valuable  because  even  approximate  appli¬ 
cation  of  the  conical-flows  method  is  tedious^  and  physical  conditions  of 
flow  in  the  tip  region  are  so  different  from  the  conditions  assumed  in  the 
theory  that  a  rough  indication  of  the  nature  of  the  theoretical  load 
should  suffice  at  the  present  time  as  a  basis  for  further  investigations. 

The  approximations  of  the  conical— flows  method^  as  it  has  been 
applied  herein,  should  be  noted:  (l)  In  canceling  lift  behind  either  wing 
half,  a  smll  amount  of  upwash  introduced  on  the  other  half  of  the  wing  - 
in  violation  of  the  boundary  conditions  —  is  ignored  and  (2),  if  the 
-trailing— edge  Mach  line  intersects  the  leading  edge,  its  reflection  will 
intersect  the  trailing  edge  and,  within  the  Mach  cone  from  the  point  of 
intersection,  only  an  estimate  of  the  loading  is  feasible.  If  the  next 
reflection  (in  the  trailing  edge)  results  in  a  second  intersection  of  the 
Mach  line  with  the  leading  edge,  the  leading— edge  correction  factor  for 
adjusting  the  two-dimensional  load  distribution  will  also  be  in  error. 
Calculations  made  to  check  the  first  of  these  approximations  showed  the 
induced-downwash  angle  to  be  less  than  half  of  1  percent  of  the  angle  of 
attack  on  the  wings  considered.  The  magnitude  of  the  inaccuracy  intro¬ 
duced  by  the  second  limitation  has  been  indicated  by  the  examples  pre¬ 
sented  herein;  the  result  is  to  limit  the  aspect— ratio  Mach  nimiber  comr- 
binations  (as  expressed  by  the  reduced  aspect  ratio  pA)  for  which  the 
conical— flows  method  will  give  accurate  theoretical  solutions.  The 
corresponding  error  in  the  leading— edge  correction  coefficient  appears, 
however,  to  be  insignificant  and  it  is  felt  that,  with  the  exception  of 
wings  tapered  to  a  point,  the  upper  range  of  reduced  aspect  ratio  PA  is 
adequately  covered  by  the  proposed  method  of  correcting  the  two- 
dimensional  load  distribution.  There  remains  only  the  case  in  which  the 
reduced  span  (P  times  the  geometric  span)  is  less  than  the  tip  chord. 

In  this  case,  the  interference  of  the  tip  flows  with  each  other  will,  in 
conjunction  with  a  subsonic  trailing  edge,  create  a  problem  not  readily 
solvable  by  either  the  conical-flows  method  or  the  two-dimensional  approx¬ 
imation.  This  problem  falls  within  the  scope  of  the  slenden-^ing  theory 
of  reference  4. 


Ames  Aeronautical  Laboratory, 

IJational  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Sept.  I6,  19^9 • 
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APEEKDIX  A 
STOBOLS 
General 

free— stream  velocity 
free— stream  Vhch  number 
i^-i 

density  of  air 

dynamic  pressure 

pressure  difference  between  upper  and  lower  surfaces,  or  local 
lift 

angle  of  attack,  radians 

circulation  on  infinite  flat  plate 

chord  of  infinite  flat  plate  (See  sketch,  p.  130 

Wing  Dimensions 

root  chord 
tip  chord 
semispan 

angle  of  sweep  of  the  leading  edge 
taper  ratio  {o^/c^) 

Coordinates 

Cartesian  coordinates  in  the  stream  direction  and  across  the 
stream  in  the  plane  of  the  wing 

coordinates  of  apex  of  oblique  trailing— edge  element 

coordinates  of  apex  of  leading-edge  element 

coordinates  of  point  on  trailing  edge,  within  tip  ^ch  cone 
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Xq  streamwise  coordinate  of  apex  of  tip  element 

Xq  largest  value  of  x^  at  which  cancellation  of  pressure  can 

affect  point  at  which  pressure  is  being  calculated 

I  streamwise  distance  of  x^y  hack  from  leading  edge^  as  a  fraction 

of  the  tip  chord  (equation  (40)) 

5^  distance  of  x^,  behind  leading-edge  tip,  as  a  fraction  of  the 

tip  chord  (equation  34)) 

Iq  distance  of  Xq  behind  leading-edge  tip,  as  a  fraction  of  the 

tip  chord  (equation  (4l)) 

Tj  spanwise  distance  of  x,y  from  leading  edge,  as  a  fraction  of 

the  chord  through  x,y  measured  perpendicular  to  the  stream 
(equation  (28)) 


In  the  following^  all  slopes  are  measured  counterclockwise  from  a 
line  extending  downstream  from  the  apex  of  the  wing  or  of  the  pertinent 
elementary  sector: 


slope  of  leading  edge  p  cot  A 
slope  of  Mach  lines 

slope  of  trailing  edge 
slope  of  ^ch  lines 

slope  of  ray  from  the  origin  ^  y 
slope  of  Mach  lines  x 

slope  of  ray  from  trailing— edge  apex 
slope  of  ^fe,ch  lines 

slope  of  ray  from  leading-edge  tip 
slope  of  ^ch  lines 


-  (Pa/m) 


slope  of  ray  from  apex  of  element  a  _  ^  J—Ja. 

slope  of  Ifech  lines  x  —  x 

a 

slope  of  ray  from  apex  of  leading-edge  element  _  ^  y—  y-^ 
slope  of  Mach  lines 

tjj  for  point  at  which  pressure  is  being  calculated 
slope  of  ray  from  apez  of  tip  element  „  y— s 

. '  . .  sr  P  . . 

slope  of  Mach  lines  x—  x„ 
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ao  the  value  of  a  corresponding  to  a  trailing— edge  element  of  which 

the  apex  lies  on  the  Mach  forecone  of  the  point  at  which  the 
load  is  being  calculated 

the  largest  value  of  a  along  which  cancellation  of  pressure  in 
the  wake  can  affect  the  point  at  which  the  load  is  being  cal¬ 
culated  through  reflection  in  the  leading  edge 

the  smallest  value  of  t^  along  which  cancellation  of  pressure 
ahead  of  the  leading  edge  can  affect  the  point  at  which  pressure 
is  being  calculated  (equation  (ll)) 

the  smallest  values  of  tg^  along  which  cancellation  of  pressure 
ahead  of  the  leading  edge  can  affect  the  point  at  which  pressure 
is  being  calculated  (equation  (16)) 


Streamwise  Components  of  Perturbation  Velocity 

basic  (uncorrected)  perturbation  velocity  as  given  by  solution 
for  triangular  wing  (equation  (l)) 

value  of  u^  at  a=0  (equation  (2)) 

constant  perturbation  velocity  on  canceling  (oblique)  sector 
in  wake 

constant  perturbation  velocity  on  canceling  sector  ahead  of  wing 
constant  perturbation  velocity  on  canceling  sector  outboard  of  tip 
symmetrical  trailing-edge  correction  to  u^  (equation  (3)) 
correction  to  u^  due  to  single  oblique  trailing*-edge  element 

(au)^  leading-edge  correction  to  u^ 

(Au-u)q  leading-edge  correction  to  Ua  due  to  reflection  of  (Au)^. 
(equation  (10)) 

(Au)g  correction  to  u^  due  to  reflection  of  leading-edge  correction 
in  wing  tip  (terms  in  E  only) 

Z^u  perturbation  velocity  induced  by  canceling  tip  effect  at  the 
trailing  edge  (equation  (47)) 
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r  .p. 


R 

C 


(AC)o 


dAC 

da 


da 


a 


a 


s 


k 

k» 

cp  or  \|r 
F(k,cp) 

K,K(k) 


E(k,(p) 

E,E(k) 

K‘(k) 


Arbitrary  Matbematical  Symbols 

real  part 

radical  term  of  leading-^dge  correction  function  (ec^uation  (9)) 

inverse— cosine  term  of  leading-edge  correction  function 
(equation  (8)) 

value  of  coefficient  of  (mx— Py)  ^  in  u/^  at  the  leading 

edge  (strength  of  the  leading-edge  singularity  in  u) 

decrement  in  the  strength  of  the  leading-edge  singularity  in 
u  due  to  reflection  of  (Au)q  at  leading  edge 

decrement  in  the  strength  of  the  leading— edge  singularity  in 
u  due  to  reflection  of  (Au)a.  at  leading  edge 

leading-edge  correction  coefficient  defined  hy  equation  (25) 

value  of  a  at  leading-edge  tip 

Elliptic  Integrals  (See  references  5  and  6.) 

modulus  of  elliptic  integral  (defined  where  used) 

the  complementary  modulus  (  </ 1  —  k^) 

argument  of  elliptic  integrals  (defined  where  iised) 

Incomplete  elliptic  integral  of  the  first  kind  of  modulus  k 
and  argument  cp 

complete  elliptic  integral  of  the  first  kind  of  modulus  k 
[K  =  F(k,|)] 

incomplete  elliptic  integral  of  the  second  kind 

complete  elliptic  integral  of  the  second  kind  of  modulus  k 
[E=E(k,|)] 

K  of  the  complementary  modulus  [K*(k)  =  K(k*)] 
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APPENDIX  B 


EVALUATION  OF  C(to)<ito  BY  SERIES  EXPANSION 

Jtq  ^■*^0 

For  integration  of  the  first  term  of  equation  (10),  the  product 
[d(Au)o/dto]  C(to)  is  integrated  by  parts  to  give 


1  d(Au)o 
dto 


C(to)dto  “  (^'^)  o 


(Bl) 


The  first  term  goes  to  zero  at  either  limit.  In  the  second  term,  t^b 
must  be  expressed  in  terms  of  x,  y,  and  to  before  taking  the  deriva¬ 
tive  of  C(to).  If  the  constant  (for  any  selected  point) 


py 

X-Co 


is  defined,  this  derivative  may  be  written 


dC(to)  _  (1+to*)  '/(To-m) (mx-py) 

d^o  (to— to*)  mco(l+m)  (l+to)  (to~^o) 


If  this  expression  and  the  expression  for  (Au)o  given  by  equation  (3) 
are  substituted  in  equation  (Bl) ,  the  result 


f  C(to)dto 

Jro 


is  obtained.  For  to  close  to  1  (small  values  of  cp),  F(  '/l-fflt^,cp)  may 

be  expanded  into  a  MacLauren  series  in  sin  cp  =  ( 1— to^ )  / (  I-hlj;^ )  as 

follows : 


uo(l+to 


K(  </l-fflt 


(Ta-ffl)(mx-3y) 

mco(l+iii) 


,r 


F(  </l-mt^,  cp) 


(to-to*)  V(l+to)( 


dti 


(B3) 
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From  the  first  tiro  terms  of  the  series,  substituted  in  equation 
(B3),  ire  may  obtain 
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APPENDIX  C 

CALCULATION  OF  TIP  LOSS  IN  LOADING  BY  CONICAL-FLOWS  METHOD 

If  the  tral ling-edge  Mach  line  intersects  the  leading  edge  on  the 
•wing,  four  components  of  the  perturbation  velocity,  as  calculated  by 
the  conical— flows  method,  must  be  canceled  outboard  of  the  tip: 3  (l)  the 

basic  velocity  u*,  (2)  the  tralling-edge  correction  (Au)^  +  /  da 


(3)  the  leading— edge  corrections  (Au)ij  =  (Aujj)^  + 


a  • 

®  d(Au) 


and  (4)  the  reflection  of  (Au)-^^  in  the  trailing  edge  (in  region 
fig,  1).  This  last  component  has  not  heen  calculated.  If  it  does  not 
extend  over  more  than  50  percent  of  the  tip  section,  the  effect  of  its 
cancellation  at  the  tip  can  prohahly  he  neglected. 

The  third  component  (Au)-j^  can  he  further  broken  down  into  the 
parts  arising  from  the  radical  term  E  in  the  leading-edge  correction 
function  and  those  arising  from  the  inverse-cosine  term  C.  The  former 
components  increase  to  infinity  at  the  leading  edge  and  will  he  treated 
in  the  same  way  as  the  basic  velocity  u^.  The  latter  remain  finite 


and  exactly  cancel  the  trailing— edge  corrections  (Au)^ 


Mu  da 
da 


at  the  leading  edge  of  the  tip  section.  Although  these  components  (the 
C  component  of  (Au)j3 ,  and  the  tralling-edge  correction)  do  not  cancel 
each  other  elsewhere  on  the  tip  section  or  in  the  stream,  the  variation 
of  each  is  fairly  small  in  the  region  affecting  the  pressures  on  the 
wing,  and,  since  the  net  Induced  effect  of  canceling  them  is  zero  and  has 
zero  slope  along  the  tip  Mach  cone,  it  will  he  assumed  that  they  cancel 
each  other  completely.  This  assumption  will  result  in  somewhat  too  low 
a  loading  near  the  trailing  edge. 

Cancellation  of  the  first  component,  the  basic  velocity  u^,  has 
heen  covered  in  reference  1.  The  result  is  a  reduction  in  u  given  by 


Au  = 


'2p(a-  y) 

x+3y 


K-2x 


'm2x2_p2y2 


F(k»,i|r) 

K' 


-  -  KE’ 


+  KE(k»,»lf) 


If  the  tip  Ilea  beyond  region  IV  (fig.  l)  further  small  components 
enter,  hut  these  will  not  he  considered  in  the  present  analysis. 
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(eqmtlon  (37)  of  reference  1),  where 
s  is  the  semispan, 

/(iii-ao){l-^ 

V  2m(aQ+l)  ^ 


J  3s(ao+m)^ 


x+P(y-s) 


K»  =  K(k«)^  and  E»  =  E(k«). 

The  part  of  the  leading-edge  correc¬ 
tion  associated  with  E(tQ)  may  he 
canceled  hy  the  following  procedure: 

For  any  one  value  of  to  (see  sketch), 
the  correction  to  u(x,y)  will  he 


dto  =  -  ^  dto  limit  rE(to,tt)  cos-1  ^•b^^c+g.^hk  + 

dt  rt2  dt  t-L,->[iL  "tc— bh 


"iL  eos-^  dti 


where  t-^  =  — Zkl  (aee  sketch) 


Integrating  hy  parts  results  in  the  elimination  of  the  first  term 
so  that 


d(Au), 


c  =  ^  ^  “o  /  ^  . 


d  t-b+tg+g^tg 


c  dt-h 

t  -t,  0 

c  h 


Equation  (C3)  can  he  integrated  analytically,  giving 


!if:ih.dt  =  2°^ 

dt  ^  jT(l+m)t^ 


(tp  —  m)  (l■^to)  (l+t-^^) 


m— t  * 

D 


where 
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To  obtain  the  total  correction  to  ■u(x,y)  for  the  cancellation  of 
the  term  in  R(tQ)  of  the  leading-edge  correction,  eq.viatlon  (04)  may  be 
integrated  graphically  by  plotting  against  (Au)^  the  remainder  of  the 
expression  in  the  right-hand  side.  An  approximate  correction  for  the 
reflection  of  the  terms  in  R(t  )  can  be  Included  by  plotting  against 
the  values  of 


Au  =  (Au)j^  + 


dAu 

da 


da 


along  the  leading  edge.  Then 


(Au)^  = 


,t-b*  = 


to=  1 


,-m)(l+to)(l+t-b*) 


dAu  (C5) 


where  and  Au  are  all  computed  for  specified  locations  of 

on  the  leading  edge  and  the  extension  of  the  leading  edge  to 
the  point  such  that  t-^*  =  —  1. 

In  practice,  equation  (C5)  appears  to  result  in  a  nearly  constant 
correction  to  the  section  loading,  as  shown  in  figures  12(a)  and  (h). 
The  work  may  he  further  simplified,  therefore,  by  performing  only  the 
relatively  simple  calculation  for  a  point  on  the  Mach  lins- 
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I  V 
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Figure  i  -  Plan  view  of  centra!  portion  of  swept- back 
wing  showing  pattern  of  Mach  tines  arising  at 
leading  and  trailing  edges. 
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Figure  2.  —  Leading-edge  eiement  and  induced-veiocity  function. 
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Figure  5.  -  Load  distribution  over  two  streamwise  sections 
of  the  tapered  swept -back  wing  calculated  by  the 
conical~flows  method . 
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L  e  a  ding  -  edge  correction. 


Oblique  tramng-edge  correction 
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25  50 
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(b)  Section  B-B;  fiy/c^  •  0.8. 

Figure  6.  -  Concluded. 


NACA  TN  1991 


o 


NACA  TN  1991 


Percent  chord 
(d)  Section  D-D-,  fiy/c^ 

Figure  7-  Conciuded. 
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Figure  8.  -  Leading-edge  velocity  coefficient  <r(x)  for  correction  of  two-  dimensional 
load  distribution. 
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Percent  chord 
(a)  Section  A- A. 


Figure  9.-  Load  distributions  on  the  tapered  wing  as  calculated 
by  the  conicah flows  method,  and  the  two-dimensional 

approximation. 
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Figure  9.  -  Concluded. 


conicahffows  method ^  and  the  two -dimensional  approximation . 


Figure  10.  -  Concluded. 


Load  distributions  on  the  untapered  wing,  m  =0.2,  as  calculated 
conical -flows  method  and  the  two-dimensional  approximation. 
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Figure  //  -  Continued. 
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Continued. 


Figure  II  -  Concluded, 
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or  94-percent  semispan. 


Figure  12.  -  Concluded. 
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